Electronic transport and vibrational modes in the smallest molecular bridge: H 2 in Pt 

nanocontacts 



Y. Garcia, 1 J. J. Palacios, 1,2 E. SanFabian, 3, 2 J. A. Verges., 4 A. J. Perez-Jimenez, 3 and E. Louis 1,2 

1 Departamento de Fisica Aplicada, Universidad de Alicante, 
San Vicente del Raspeig, Alicante 03690, Spain. 
2 Unidad Asociada UA-CSIC, Universidad de Alicante, 
San Vicente del Raspeig, Alicante 03690, Spain. 
s Departamento de Quimica-Fisica, Universidad de Alicante, 
San Vicente del Raspeig, Alicante 03690, Spain. 
* Departamento de Teoria de la Materia Condensada, 
Instituto de Ciencias de Materiales de Madrid ICMM-CSIC, Cantoblanco, Madrid 28049, Spain. 

(Dated: February 2, 2008) 

We present a state-of-the-art first-principles analysis of electronic transport in a Pt nanocontact 
in the presence of H2 which has been recently reported by Smit et al. in Nature 419, 906 (2002). 
Our results indicate that at the last stages of the breaking of the Pt nanocontact two basic forms of 
bridge involving H can appear. Our claim is, in contrast to Smit et al.'s, that the main conductance 
histogram peak at G s=s 2e 2 /h is not due to molecular H2, but to a complex Pt2H2 where the 
H2 molecule dissociates. A first-principles vibrational analysis that compares favorably with the 
experimental one also supports our claim. 

PACS numbers: 73.63.Fg, 71. 15. Mb 



Prediction ability is, possibly, the most important at- 
tribute associated with first-principles numerical imple- 
mentations of quantum transport theory in molecular 
and atomic-scale systems. This desirable property re- 
lies on the correct handling of two closely-related factors: 
The atomic structure and the electronic structure of the 
nanoconstriction or molecular junction. The calculation 
of the electronic structure in current-carrying situations 
is a topic of increasing interest and has received due at- 
tention in recent years [J S H S H H 

The compu- 
tation of the atomic structure, which is well-established 
on standard geometry relaxation methods, is, however, 
much more demanding from a first-principles computa- 
tional point of view and it is not subject to scrutinity on 
a regular basis in molecular electronics. Besides, there 
are inherent uncertanties associated with the various ex- 
perimental procedures that can only be addressed by sta- 
tistical analysis. 

In this paper we show how a combination of state-of- 
the-art current-carrying first-principles electronic struc- 
ture calculations and a careful structural analysis can 
shed light on a recent experiment :8J,j9| of electronic trans- 
port in Pt nanocontacts where, apparently, H2 molecules 
anchor themselves to the metal and modify the standard 
conductance histogram of Pt. A H2 molecule is the sim- 
plest molecule and is ideal to carry out a detailed struc- 
tural, vibrational, and transport analysis. On the basis 
of the extensive ab initio calculations we have performed, 
our main conclusion is the following: Contrary to the au- 
thors' explanation for the conductance histograms where 
the main peak around G ~ Go (for Go = 2e 2 /h) is at- 
tributed to a configuration where a H 2 molecule bridges 
the gap between Pt electrodes (see also Ref. 0) , we at- 



tribute this peak to a situation where the H2 molecule 
dissociates at the nanocontact and forms a stable com- 
plex Pt2H-2 where the two H atoms chemisorb between 
electrode tip atoms [see left inset in Fig. Our pro- 
posed structure presents vibrational modes in accordance 
with the reported inelastic effects and with a more re- 
cent vibrational analysis 9J . Finally, we speculate that 
the bridge with H in molecular form could account for 
the peak around G ~ 0.1 Go reported in Ref. 0, but not 
commented on there. 

For the structural and vibrational analysis of the 
Pt nanocontact with H we have used the GAUS- 
SIAN98 code[m and for the transport calculations we 
have employed the Gaussian Embedded Cluster Method 
(GECM), recently developed by the authored S E3 
which, in turn, makes use of GAUSSIAN98. The GECM 
consists of a first-principles evaluation at the density 
functional theory (DFT) level of the density matrix and 
the Fock matrix F (or sclfconsistent hamiltonian) of the 
region determining the transport (here the H atoms and 
the Pt electrode atoms in proximity to them) in a current- 
carrying situation. As regards the DFT calculations for 
the core clusters (see insets in Figs. Q] and |3J| we used the 
B3LYP ninctional(n| , th e basis sets and core pseudoten- 
tials described in Ref. El El El for Pt and the cc-pVTZ 
basis for rlEU • Here we are concerned with the linear 
response conductance which is calculated in the standard 
manner 1 16|: 

G = ^Tr[th}. (1) 
In this expression, Tr denotes the trace over all the or- 
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bitals in the cluster. The matrix t is given by 



(2) 



where the Green function matrix (in orthogonal basis) is 



G (±) = \(E±iS)l -F- S (±) 



(3) 



while 1 is the unity matrix. 
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The matrices and r L are 
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respectively, 
JR -r ^ L represents the 
contribution of the semiinfinite left and right leads not 
included in the self-consistent evaluation^ In order 
to single out the contribution of individual channels to 
the current we diagonalize the matrix t't. 

As a starting point we consider the conductance of 
a clean Pt nanocontact without H 2 . It would be de- 
sirable to simulate the whole breaking(formation) pro- 
cess of the nanocontact on stretching(indentation). This 
is, unfortunately, computationally not feasible from first 
principles for Pt. Out of the many possible atomic ar- 
rangements, we have chosen to study the conductance 
for the following nanocontact atomic configurations: 9- 
4-1-4-9 (single-atom contact) and 9-4-1-1-4-9 (chain con- 
tact) for the (001) direction and 6-3-1-3-6 (single-atom 
contact) and 6-3-1-1-3-6 (chain contact) for the (111) di- 
rection (the numbers indicate the number of atoms per 
plane). We have considered two possible plane order- 
ings for the (111) direction: 6B-3A-. . . -3A-6B and 6C- 
3A-. . .-3C-6A (remember that the bulk one is . . . A-B- 
C-A-B-C . . . ). These nanocontacts correspond to cases 
with and without mirror symmetry with respect to the 
middle plane perpendicular to the stretching axis z [the 
(001) nanocontact has also been chosen with mirror sym- 
metry] . Specifically, we have taken bulk atomic distances 
for the atoms in the outer planes, but we have performed 
full relaxation in all the other coordinates, including the 
distance between the two outer planes. The resulting 
structures are thus representative of the ones respon- 
sible for the last (first) conductance plateau on stretch- 
ing(indentation) cycles. 

The conductance reaches ~ 2Q at the Fermi energy 
(set to zero) for the two short-chain cases with mirror 
symmetry and drops slightly for the other one [see Fig. 
ffl^b)]. Two conduction channels mainly contribute to 
the conductance in the three clusters considered with 
the two-atom chain [see Fig. Efd)]. Given the polycrys- 
tallinc nature of the electrodes, one should not expect 
mirror symmetry to be the dominant situation in the 
experiment so we can expect an average conductance 
around Q w 1.75Ga- This value nicely coincides with 
the position of the lowest conductance histogram peak 
reported in Ref. 0- The conductance of the single-atom 
contacts presents similar values for the (111) direction 
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FIG. 1: (a) Conductance as a function of energy (referred 
to the Fermi energy) for the (lll)-oriented 6-3-1-3-6 clusters 
with (blue line) and without (red line) mirror symmetry and 
for the (lOO)-oriented 9-4-1-4-9 cluster; (b) same as in (a) but 
for the (lll)-oriented 6-3-1-1-3-6 clusters with (blue line) and 
without (red line) mirror symmetry and for the (lOO)-oriented 
9-4-1-1-4-9 cluster; (c) channel transmission for the 6-3-1-3-6 
cluster with mirror symmetry [see inset in (a)] ; (d) same as in 
(c) but for the 6-3-1-1-3-6 cluster [see inset in (b)] with mirror 
symmetry. 



(although, as Fig. Ofc) shows, three channels now con- 
tribute to transport), but is appreciable higher in the 
(001) direction [see Fig. ma)]. Although we cannot af- 
ford a statistical analysis, the fact that the position of 
the lowest conductance histogram peak shifts to higher 
values when single-atom contacts are favoured (on only 
considering indentation cycles or higher bias voltages as 
reported in Ref. Il7|) is consistent with our results. The 
conductance is also here noticeably dependent on the de- 
tailed atomic arrangement of the electrodes (note that 
the conductance increases appreciably on the scale of 1 
eV below and above the Fermi energy). This is due to the 
fact that the Fermi level lies on the edge of the d band 
for bulk and the contact atom has almost bulk coordina- 
tion. This is also fairly consistent with the fact that Pt 
histo gram peaks are broad in comparison to, e.g., those 
of Au[3- 

When the same conductance measurements are per- 
formed in a H2 atmosphere the histograms change no- 
ticeably. Inelastic vibrational spectroscopy along with 
theoretical support lead the authors in Ref. |g to conclude 
that a H2 molecule can position itself along the stretching 
axis of the Pt nanocontact right before it breaks apart. 
The authors claim that this configuration accounts for 
the main histogram peak for conductance where G ~ Go- 
With this claim in mind we consider nanocontacts like 
that shown in the inset of Fig. IHb), but with a H2 
molecule anchored in between electrode tip atoms. The 
results presented below are for the smallest cluster in 
the (001) direction, 4-1-H2-1-4, but we have performed 
similar calculations for (111) nanocontacts for reassure- 
ment. We maintain bulk interatomic distances for the 
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four atoms in the two outer planes whose variable respec- 
tive distance d mimics the effect of the piezo voltage. For 
every d we let the positioning of the H atoms and that of 
the two contacting Pt atoms relax, starting from a con- 
figuration where the H2 molecule sits along the stretching 
axis. Full blue circles in Fig. |21 represent the total energy 
vs. d for all the stable configurations we have found in 
which the molecule stays on the z-axis (see right inset in 
Fig.GJ. Beyond d 9.5A the molecule sticks to one of the 
electrodes or dissociates. Below d « 8.25A the molecule 
tends to loose its alignment with the z-axis. (As a guide 
to the eye we also present with empty circles the energy 
for configurations where the molecule is forced to stay 
along the z-axis). In all these points the H2 retains its 
molecular form. At the middle point (d ~ 9 A) the H-H 
bond distance and the H-Pt distance turn out to be sim- 
ilar to those proposed in Ref. |g. Figure^a) shows Q(E) 
for d = 9.25A ( we chose the distance that more closely 
reproduced the vibrational mode frequency reported in 
Ref.0, see Table I). The conductance at the Fermi energy 
turns out to be remarkably smaller than Go- A similar 
calculation for the whole range of stable configurations 
shown in Fig.|5]gave a conductance in the range 0.2-0.5C?o- 
To our ease this result agrees with the intuitive picture: 
Closed-shell molecules cannot conduct if the molecular 
character is maintained after contact with the electrodes 
and the charge transfer is much smaller than one. Both 
conditions apply here. The charge transfer from the H2 
molecule to the electrodes is typically w 0.15 or smaller 
and the bonding and antibonding molecular levels main- 
tain their character since they are clearly visible in the 
density of states (not shown here) around -6 eV and 17 
eV, respectively. A value of conductance close to Go is 
still possible through hybridization of the Pt states with 
those of the molecule In our calculations, however, 
this does not occur at the Fermi energy, but at E « —1.5 
eV. We believe that this H2 bridge actually accounts for 
a secondary peak around Q « 0.1C/o not discussed, but 
clearly present in the reported dataQ- 

So, what accounts for the main conductance peak at 
Qo?. There is little doubt that this peak is H-related, 
but here we propose an alternative to the H2 molecu- 
lar bridge discussed above. We now perform calcula- 
tions similar to the previous ones, but initially placing 
the molecule axis perpendicular to the stretching z-axis. 
We let the Pt tip atoms and the H atoms relax towards 
equilibrium for several values of d without assuming the 
H2 to be in molecular form. The total energy of the sta- 
ble configurations we have found are represented by full 
black circles in Fig. |2 (Again, for completeness, we also 
present results when the H atoms are not allowed to leave 
the plane perpendicular to the z-axis). The two H atoms 
are only in molecular state at the upper limit of this 
curve (d « 8 A). For smaller values of d the H2 molecule 
breaks apart and the H atoms position themselves up to 
2.0 A apart from each other forming a double atomic H 
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FIG. 2: Total energy of a H2 molecule in a Pt nanocontact 
versus distance between electrodes (outer planes). Blue cir- 
cles correspond to the H atoms positioned along the stretching 
axis (z) in molecular form (see right inset) and black circles 
correspond to the H atoms positioned on the plane perpen- 
dicular to the stretching axis (see left inset). Empty symbols 
have been added as a guide to the eye which correspond to 
restricted relaxations where the H atoms are forced to stay on 
the z-axis (blue) or on the plane perpendicular to it (black). 
Arrows indicate configurations for which vibrational mode fre- 
quencies closest to the experiment are obtained. 

bridge as the one shown in the left inset of Fig. The 
results of conductance for d = 6. 5 A (a distance for which 
vibrational modes are closer to the experimental ones, 
see below and Table I) are depicted in Fig. IHIb). The 
conductance at the Fermi energy approaches closely Go- 
The eigenchannel transmission results are also consistent 
with the noise analysis in Ref. since, the fact that fluc- 
tuations do not strictly vanish at Go m ay indicate that 
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FIG. 3: (a) Conductance as a function of energy (referred to 
the Fermi energy) for the cluster shown in the inset where a H2 
molecule sits along the stretching axis of the Pt nanocontact 
(d = 9.25A). Dotted lines in (c) represent the transmission of 
the individual eigenchannels for the cluster in (a) . (b) Same as 
in (a) but for the structure shown in the inset where the two H 
atoms form a double atomic bridge (d = 6.5A). Dotted lines in 
(d) represent the transmission of the individual eigenchannels 
for the cluster in (b). The distance d was chosen to give, in 
both cases, the main vibrational mode frequency observed in 
the experiments. 
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more than one channel is contributing to the conduc- 
tance. At this point the distance between tip Pt atoms 
is approximately that of bulk Pt and that between tip 
Pt atoms in Fig.QIb). When comparing with the results 
presented in Fig. ^b) we see that the H seems to block 
out one of the channels. It is important to finally stress 
here that the exact value of the conductance is not very 
dependent on the basis set and the functional. However, 
the atomic configuration of the electrodes in the cluster 
and the model used to represent the bulk of these elec- 
trodes introduces a source of uncertainty in the results. 
We believe that this is the source of the discrepancy be- 
tween our results and those presented in Refs. @ and lift 
After considering various other possibilities in terms of 
cluster sizes and orientations, we can safely conclude that 
the stark difference between the conductance of the two 
basic bridge structures studied in this work is robust. 

We have also performed a full analysis of the H-H vi- 
brational modes in a configuration in which the H2 axis 
is perpendicular to the z-axis and for the distance that 
gives a conductance close to one quantum (d = 6. 5 A see 
Figs. [5] and |3J. The results for the mode frequencies are 
reported in Table I. The first mode shown in the table (H 
atoms moving in opposite directions along the stretching 
axis) is the one that has a frequency closer to the exper- 
imental data (a frequency of 516±21 cm -1 was reported 
in 0). The isotopic effect shows up clearly with frequen- 
cies 439 and 386 cm -1 for H-D and D-D, to be compared 
with those obtained by multiplying the H-H frequency 
by the mass ratii (445 and 392 cm" 1 , respectively). It 
is worth mentioning that, although for D-D the isotopic 
effects shows up for all modes, this is no longer true for 
H-D. This mode induces a large variation of the conduc- 
tance (see Table I), again, an effect that is not the same 
for all modes. Note also that the mode with the two H 
atoms moving in the same direction along the stretch- 
ing axis, which also induces a significant change in the 
conductance (see second row in Table I), has a frequency 
close to that of a second mode observed in the experi- 
ments (around 1200 cm" 1 , see Ref. 0). The latter is not 
reproduced by configurations with a molecular H2 bridge 
(see Table I). 

In conclusion we have analyzed the conductance and 
vibrational modes of a Pt nanocontact in the presence 
and absence of H. Our results for clean Pt agree with 
most data reported in the literature. On the contrary, our 
results on Pt with H question the current interpretation 
of the experiments reported in Ref. @. While H 2 seems 
to block the current, a novel complex Pt2H2 formed at 
the neck could be responsible for the reported histogram 
peak near Go . The vibrational modes of this complex also 
agree with the experimental ones. 
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and No. MAT2002-04429-C03 and from the Universidad 
de Alicante. We thank C. Untiedt for sharing with us his 



TABLE I: Frequencies (in cm" 1 ) of the six H-H vibra- 
tional modes in the configuration for d = 6.5A where the 
H2 molecule dissociates and places itself with its axis perpen- 
dicular to the stretching axis (see Fig. |2J . Arrows indicate 
atomic displacements in the plane defined by H and tip Pt 
atoms while ± denotes displacements perpendicular to this 
plane. In the last column the numerical results for the abso- 
lute value of the relative change in conductance induced by 
a 0.1 A displacement of the H atoms are reported. The H-H 
vibrational modes for the configuration with the H2 molecule 
placed along the stretching axis and d = 9.25A are shown in 
parenthesis. 
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341 (1707) 


307 
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0.012% 



point of view on the experimental results. 
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